Weak measurement of the Dirac distribution 
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Recent work [J.S. Lundeen et al. Nature, 474, 188 (2011)] directly measured the wavefunction by 
weakly measuring a variable followed by a normal (i.e. 'strong') measurement of the complementary 
variable. We generalize this method to mixed states by considering the weak measurement of the 
product of the two observables, which, as a non-Hermitian operator, is normally unobservable. 
This generalized method provides mixed states an operational definition related to the operator 
representation proposed by Dirac. Uniquely, it can be performed 'in situ', determing the quantum 
state without destroying it. 
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The wavefunction ^ is at the heart of quantum me- 
chanics, yet its nature has been debated since its incep- 
tion. It is typically relegated to being a calculational 
device for predicting measurement outcomes. Recently, 
Lundeen et al. proposed P, a simple and general oper- 
ational definition of the wavefunction based on a method 
for its direct measurement: "it is the average result of a 
weak measurement of a variable followed by a strong mea- 
surement of the complementary variable." The 'wave- 
function' referred to here was introduced along with the 
Schrdinger Equation and can be thought of as a special 
case of a quantum state, known as a 'pure state.' The 
general case is represented by the density operator p, 
which can describe both pure and 'mixed' states. The 
latter incorporates both the effects of classical random- 
ness (e.g., noise) and entanglement with other systems 
(e.g., decoherence) . The density operator plays an im- 
portant role quantum statistics, quantum information, 
and the study of decoherence. As quantified by the Pu- 
rity fi — Tr[p^] < 1, a density operator can range from 
pure (/i = 1) to completely mixed (/i — 1/N, where N is 
the Hilbert space dimension). We investigate what be- 
comes of our operational definition of the wavefunction 
when it is applied to general quantum states. 

The standard method for experimentally determining 
the density operator is Quantum State Tomography Q. 
In it, one makes a diverse set of measurements on an 
ensemble of identical systems and then determines the 
quantum state that is most compatible with the measure- 
ment results. An alternative is our direct measurement 
method, which may have advantages over tomography, 
such as simplicity, versatility, and directness. As com- 
pared to tomography, which works with mixed states, a 
significant limitation of the direct measurement method 
is that it has only been shown to work with pure states. 
The most common goal of state determination is to eval- 
uate a system's difference from a target quantum state, 
which might be, for example, a potential resource for a 
quantum information protocol. Since the target is usu- 
ally a pure state, this difference often due to mixedness, 
and so it can be quantified by the purity of the state. 

In this paper, we consider whether our direct method 



can be generalized to density operators. Previous 
works have developed direct methods to measure quasi- 
probability distributions such as the Wigner function [4] , 
Husimi Q-function and the Glauber-Sudarshan P- 
function [6[. These are position- momentum (i.e. 'phase- 
space') distributions that are equivalent to the density 
operator, and have many, but not all, of the proper- 
ties of a standard probability distribution. The most 
used quasi-probability distribution is the Wigner func- 
tion Wp{x,p), which has the following properties: (1) it 
is real; (2) less than one, \Wp\ < 1; (3) its marginals give 
correct predictions for position x and momentum p prob- 
ability distributions (i.e. J Wp{x,p)dx = Prob(p)); and 
(4) the mean value of an observable A is just an overlap, 
(A) = 2iTh J J Wp{x,p)WA{x,p)dxdp. However, unlike a 
standard probability Wp can be negative. This is what 
enables it to be compatible with the predictions of quan- 
tum mechanics. The Wigner function can be directly 
measured Q by displacing the system in phase space and 
then measuring the parity operator (this is a nontrivial 
requirement, see Q). Equivalent ly, the integral of the in- 
terference between a pair of rotated and displaced repli- 
cas of the system will give the Wigner function The 
Husimi Q-function can be directly measured by an eight- 
port homodyne apparatus or by projection on the har- 
monic oscillator ground state [lOj . These phase-space dis- 
tributions are created to be the closest quantum analogs 
to a classical probability distribution. In this sense, they 
are inherently amenable to direct measurement. 

A lesser known phase-space distribution that can be 
used to represent a quantum operator was introduced by 
Dirac in 1945 in his paper, "On the Analogy Between 
Classical and Quantum Mechanics." The Dirac distribu- 
tion has since been extended to discrete Hilbert spaces 
11| . In various guises it has been investigated periodi- 



cally during last half-century ^U]. In optics, variations of 
the Dirac distribution have been used widely, appearing 
in Walther's definition of the radiance function in radiom- 
etry [l^ and Wolf's specific intensity (as pointed out 
in Ref. [ll|). Dirac showed that an operator A could be 
represented by its overlap with the basis states of two 
non-commuting variables. Specifically, one can represent 
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an operator in phase-space as Sa{x,p) = {x\ A \p) ■ {p\x). 
li A = p, the Dirac distribution is representation of the 
quantum state of a system and shares features (2) to 



(4), from above, with the Wigner function Unhke 



the Wigner function, it is compatible with Bayes' law 
and, thus, is consistent with a quantum analog of classi- 
cal determinism [l5| . Hence, it has many of the desired 
features of quasi-probability distributions but also has a 
novel logical consistency. 

However, the Dirac distribution has two peculiar fea- 
tures: one, it is sensitive to the ordering of the non- 
commuting variables (i.e. (a;|A|p) • {p\x) ^ {p\A\x) ■ 
{x\p))] and two, Dirac noted that although his distribu- 
tion "was developed to provide a formal probability" for 
X and p "it turns out to be in general a complex number." 
Surprisingly, Dirac saw the sensitivity of his distribution 
to the arbitrary choice of ordering not as a drawback 
but, rather, as a desirable characteristic, since it em- 
phasized the key difference between quantum and classi- 
cal mechanics: non-commutivity of operators [l6| . Even 
accepting this, the Dirac distribution's complex nature 
has likely inhibited its acceptance as a quasi-probability 
distribution. Other complex phase-space operator rep- 
resentations have proven useful in quantum optics [17| . 
And ways in which complex probabilities can be made 
logically consistent with a frequentist interpretation of 
probability have been outlined [ISj . Surprisingly, despite 
being complex we will show that the Dirac distribution 
is directly measurable, much like the Wigner function. 
Notably, in our method, the choice of variable ordering 
has a clear physical meaning. Thus, it removes a degree 
of arbitrariness from the distribution's definition. 

We begin by considering what happens to our method 
for directly measuring the wavefunction when the state 
is not pure. At the heart the direct method is weak 
measurement, which we now review. Upon a standard 
(i.e. 'strong') measurement of a quantum system, the den- 
sity operator changes depending on the measurement 
outcome. This change can be viewed in two ways: as 
a straightforward update of our knowledge of the state 
of the system or, at least in some situations 18|, as a 
physical disturbance due to the interaction between the 
system and the measurement apparatus (e.g. the Heisen- 
berg Microscope [Tijl). To perform a weak measurement 
one reduces the strength of the interaction between the 
measured system and the measurement apparatus. This 
results in a commensurate reduction in the amount of in- 
formation one receives from the measurement and, also, 
a reduction in the disturbance induced by the measure- 
ment. This tradeoff is inherent in quantum mechan- 
ics: measurement precision for system disturbance pot . 
While a weak measurement on a single system provides 
little information, by repeating it on an arbitrarily large 
ensemble of identical systems one can precisely determine 
the average measurement result with arbitrary precision. 
This average is simply the standard quantum expectation 



value, (^fj A l^f), where A is the observable measured and 
1 4') is the pure state of the system. This is true indepen- 
dent of the strength of the measurement. 

Over the last decade, interest in weak measurement 
has grown as researchers have realized its potential for in- 
terrogating quantum systems in a coherent manner 21 1 . 
Weak measurement theory has been used to model and 
understand photonic phenomena in birefrigent photonic 
crystals [ll], fiber networks [131, cavity QED and 
quantum tunelling |25|. Weak measurement provides in- 
sight into a number of fundamental quantum effects, in- 
cluding the role of the uncertainty principle in the double- 



slit experiment 3, 2^, the Legget-Garg inequality 27 1, 
the quantum box problem [il], and Hardy's paradox 29 1. 
Weak measurement has also been used to amplify small 
experimental effects fsO^ and as feedback for control of 
a quantum system [31|]. Weak measurements have been 
demonstrated in both classical [12] and non-classical sys- 
tems 



A distinguishing feature of weak measurement is that, 
in the limit of zero interaction, the quantum state of the 
system remains unchanged. Subsequent measurements 
can now provide additional information about that ini- 
tial quantum state. Consider a subsequent strong mea- 
surement of observable C that results in outcome c (cor- 
responding to eigenstate eigenstate |c)). The average re- 
sult of the weak measurement of A in the sub-ensemble 
of systems giving C = c is called the 'Weak Value' and 
is given by [34j] . 



{c\A\^) 

(cm 



(1) 



Surprisingly, the weak value can be outside the range o f 
the eigenvalues of A and can even be complex 35-33- 
The concept of weak measurement is universal |38l440| | 
but insight into the meaning of complex weak values can 
be gained by considering a p articular model of measure- 
ment by Von Neumann |4l| . In this general model, the 
measurement apparatus has a pointer that shifts in po- 
sition to indicate the result of a strong measurement of 
A. In weak measurement, not only does the average po- 
sition of pointer shift but also the average momentum 
of the pointer. These shifts are proportional to the real 
and imaginary parts of the weak value, respectively, thus 
giving them straightforward physical manifestations 35- 
37[. The complex nature of the weak value is what en- 
ables us to directly measure real and imaginary parts of 
the wavefunction and, as we show later, directly measure 
the Dirac distribution. 

We now review our method for the direct measurement 
of the wavefunction. The concept is general, however 
here we consider the case of a discrete Hilbert space. In 
this space, one is free to choose the basis {|a)} (associated 
with observable A) in which the wavefunction will be 
measured. The method consists of weakly measuring a 
projector in this basis tTq = \a) (a|, and post-selecting on 
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a particular value 60 of the complementary observable B. 
By 'complementary' we mean that (a|6o) — l/V^ for all 
a, where A'' is the dimension of the Hilbert space. That 
is, the overlap is real and constant as function of a. The 
existence of state |6o) is guaranteed by the existence of at 
least two mutually unbiased bases (MUB) in any Hilbert 
space 
by 



42l |. Using Eq. ([T]), the quantum state Y^) is given 



If) = V 



(2) 



where (tTq)^ is the weak value and w is a constant that is 
independent of a. Thus by stepping through the values 
of a in a series of weak measurements one can directly 
measure I^E") represented in the a basis. 

To extend our method to mixed states we need to cal- 
culate the weak value of a system described by a density 
operator. There always exists an orthonormal basis {|A)} 
in which, 



(3) 



Thus, a mixed state can always be thought of the random 
preparation of a set of orthogonal pure states |A), each 
with probability p\. Refs. [43| showed that the weak 
value will simply be a weighted sum of the weak values 
for every pure state in the decomposition given in Eq. 



(4) 



where P (A|c) is the probability of the initial state being 
prepared in |A) given the system was post-selected in |c). 
Using B ayes' Theorem to simplify this (as in [i^) we 
arrive at: 



{c\Ap\c) 
{c\p\c) ■ 



(5) 



And applying this to our direct measurement method we 
find. 
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Examination of Eq. (jB)) shows that only 27V real param- 
eters are found by scanning a. This will not generally 
be sufficient to determine all the parameters in p, which 
has iV^ — 1 real parameters. As might be suspected, our 
method for the direct measurement of the wavefunction 
cannot be used to determine a mixed state. 

We now consider whether the direct measurement tech- 
nique can be used for a more modest task: to deter- 
mine whether a state is mixed. This could be used as 
an indicator for when the method should not be used. 
Or, if it is also a measure of state purity ^, it could 



be used to determine the quality of prepared quantum 
states. The weak value in Eq. (jG)) is equal to sum of 
one row of the density matrix, Paia2 — {'^i\p\o-2)i i-e. 
^'^"■i)^P = N-Pbaba Paias, where Pbofoois the probability 
that the measurement oi B — bo. In effect, one measures 
a distinct weak value for each pure state in the statistical 
decomposition in Eq. This simply adds noise to our 
determination of the weak value for the density operator. 
In the limit of zero interaction strength, the weak mea- 
surement exhibits maximal statistical uncertainty and so 
additional noise would be imperceptible. Thus, we con- 
clude that the original direct measurement method does 
not contain any signature of the purity of the state. 

We now consider what happens if one replaces the 
strong measurement of B with a weak measurement. 
Specifically, we investigate the weak measurement of 
the product of projectors from the two MUB, Sab = 
|6)(6|a^^a|. Since there is no longer any post-selection, 
this average is calculated in the same manner as the stan- 
dard quantum expectation value: 



{Sab)p = Tr [Sabp] 

= {a\p\b) {b\a} 
= Sp{a,b), 



(7) 
(8) 
(9) 



where Sp{a, b) is the discrete Hilbert space version of the 
Dirac distribution as defined in 11] (there exist anal- 
ogous discrete Wigner functions 44]). This is the main 
result of the paper; A joint weak measurement of comple- 
mentary variables A and B directly measures the Dirac 
distribution at values a and b. By scanning a and b, so 
as to measure the Dirac distribution for over all values 
of (a, 6), one completely determines the density operator. 
Similarly, a joint weak measurement of Sxp = \p) (p|a;) (a;| 
on a mixed state p gives the phase-space version of the 
Dirac distribution, Sp{a,b). 

In order to calculate the density operator from 
the Dirac distribution one must know {b\a) = 
exp{i9ab) /VN- Unfortunately, other than in a select 
few cases, it is not generally known what are the bases in 
the MUB set (for any given Hilbert space). Thus, a gen- 
eral formula for 9ab is also unknown. However, if {|a)} is 
taken to be the standard basis (i.e. J2a=o (^1 ~ ^' 
the identity operator) then one MUB, which we take 



to be {|6)}, will always be the Fourier basis 

=0 Ifl) 6xp (i27ra6/iV) /-s/TV. In this case, 9ab — 
—2'iTab/N , where a and b are integers solely used to enu- 
merate the states such that 0<a, & <A^— 1. With these 
choices for our complementary bases the density operator 
is simply related to the Dirac distribution by a Discrete 
Fourier Transform, 



N-l 



Oa,a, = Sp{ai,b)e'^^'^^^- 



■02)/Af 



(10) 



b=0 



This explicitly shows that average result of the joint weak 
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measurement, {Sab)p, contains the same information as 
the density operator. 

However, there is no need to ever transform to the 
density operator as all of the system's properties can be 
found directly with the Dirac distribution, Sp{a,b). It 
can be used to directly calculate the average result of 
measuring an observable O on state p: 



N-l 

{0)^N-Y, Sp{a,b)So{a,b)* 

a,b=0 



(11) 



where So (a, b) can be calculated by Eq. d?]), and * is the 
complex conjugate. With O chosen to be the identity 
operator /, we find that X]^6=o ^) ~ [p] = 1- 

In other words, the Dirac distribution is normalized in 
the same manner as a probability distribution. With 
O = p, we find the purity p = N ■ J2aM=o («' ^)l^ ' 
which reaffirms that purity is a global property of the 
density operator and thus, we are unable to measure pu- 
rity without completely determing p. And finally, with 
O = \a) {a\ or \b) {b\, we see that the a and b marginals 
are equal to the probability distributions of outcomes a 
and b, e.g. J2a=o ^) = Prob(6). Consequently, the 

result of our joint weak measurement, the Dirac distri- 
bution of the density operator, is a capable alternative 
to the standard quantum quasi-probability distributions, 
such as the Wigner function. Its peculiarity is that it is 
complex whereas probabilities are real. Nonetheless, our 
method for directly measuring it provides an operational 
meaning to both its real and imaginary parts; they ap- 
pear right on our measurement apparatus, in the shifts 
in the two conjugate observables of the pointer, e.g., x 
and p. 

One cause for concern in our method is that Sab is not 
Hermitian (i.e. = F, where f indicates the Hermi- 
tian Adjoint). According to the postulates of quantum 
mechanics this means it can not be an observable [46| . 
The typical justification for this postulate is that only a 
Hermitian operator has real eigenvalues. This is critical 
for two reasons: One, Hamiltonians containing these ob- 
servables must be Hermitian in order to conserve prob- 
ability. This argument has been challenged recently in 
regards to PT symmetry, And two, because we only 
ever record real values in experiments. An established 
counter-example to this justification is the weak value, 
Eq. ([1]) , which can be complex and yet can be measured 
straightforwardly. One way to view the familiar com- 
plex values in quantum mechanics (e.g., in the Schrdinger 
equation) is as a mechanism for mathematical simplifica- 
tion, enabling the combination of coupled equations into 
one equation. In the spirit of this view, the complex weak 
value is a just a simplified expression for two real aver- 
age measurement outcomes. Complex eigenvalues create 
no inconsistencies in this simplified formulation. Con- 
sequently, instead of Hermitian we contend that weakly 
measured observables are only required to be Normal, 



The above discussion assumes that the weak measure- 
ment of the non-Hermitian operator Sab is actually pos- 
sible. It is reasonable to question whether, in principle, 
there is a measurement apparatus that can jointly im- 
plement a weak measurement of \a) {a\ and \b) The 
answer is yes. One possibility is to have the same mea- 
surement pointer coupled to both A and B variables [3] . 
Another possibility is to weakly measure a and b projec- 
tors separately, coupling each to its own measurement 
pointer, and then evaluate the average correlations in 
the two pointers (e.g., {xiP2) and {X1X2)) H, 49 1. Both 
options have been demonstrated experimentally 2^, l50| . 
Either possibility allows one to find the average result of 
the weak measurement of Sab- 

We now return to the operator ordering in the Dirac 
distribution, which as it stands is chosen arbitrarily. An 
examination of Eq. ([7]) reveals that in our definition of 
Sab the measurement of A precedes the measurement of 
B. One might hypothesize that the order in which two 
weak measurements conducted would not matter since 
the wavefunction is nominally unchanged by both mea- 
surements. In Ref. [4^, the authors considered sequen- 
tial weak measurements of observables E and F and three 
possibilities were distinguished: measure E then F; mea- 
sure F then E; and, measure E and F simultaneously. 
By analyzing each of these with the Von Neumann model 
of measurement, they found that each had a different 
weak value if E and F do not commute. In the con- 
text of our direct measurement of the Dirac distribution, 
the operator ordering choice is, thus, set by the choice of 
time-ordering in which one conducts the measurements 
of A and B. 

Experimentally implementing a joint weak measure- 
ment can be challenging since the signal to noise of 
the weak measurement of n operators SNR(n) goes as 
SNR(l)". We now outline an easier method to deter- 
mine the density matrix. It consists of scanning the 
post-selected state in our original method for the direct 
measurement of the wavefunction. In this scenario, we 
replace 60 with b in Eq. ([5]), where {\b)} is the Fourier 
basis. The weak value is given by, P (6) (Tra)^ = Sp{a, 6), 
where the probability to post-select B = b is P (b) = ptb- 
To consider the experimental meaning of the left-hand 
side of this equation, we return to the general formula 
for the weak value, Eq. ([5]). In the von Neuman model, 
the pointer shift, which indicates the weak value of a 
general observable A, is determined by measuring the 
average of some pointer observable D = J^i'^i Mi) (dil 
(e.g. D = X) in the subset of systems with outcome 
C = c or = kJ2idiP {di\c), where fc is a pro- 

portionality constant accounting for the strength of the 
measurement. Then using Bayes' Theorem, we find 
(c|Ap|c) — kJ2idiP {di,c), where P{h,g) means the 
probability of getting both h and g outcomes in one trial. 
Applying this result to the weak measurement of \a) {a\ 
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while post-selecting B = b we find, 

Sp{a, b) = k'^^diP {dt,c) 



(12) 



Consequently, again we can directly measure the Dirac 
distribution, though not through the usual evaluation of 
the weak value. Instead one records the joint probabil- 
ity for outcome c on the system of interest and outcome 
di on the pointer, when weak measuring \a) {a\. In this 
way we can completely determine the density operator 
through Eq. (fTO|). This determination has a key advan- 
tage over standard quantum state tomography in that it 
only requires measurements in two of the system's bases. 

In summary, by measuring a pair of complementary 
variables of a system it is possible not only to directly 
measure its wavefunction but, also, to determine its den- 
sity operator. To measure over the extent of wavefunc- 
tion we only need to scan the first variable. To mea- 
sure the density operator we must, additionally, scan the 
second variable. Whereas the first variable must be mea- 
sured weakly, the second can either be strongly or weakly 
measured. The latter option corresponds to the measure- 
ment of the product of projectors for two complementary 
observables, such as x and p. This product operator is a 
non-Hermitian operator and thus has no corresponding 
observable relevant for standard measurement. Nonethe- 
less, it is weakly measurable, the average result of which 
is a two-dimensional complex distribution that was first 
introduced by Dirac. With the distribution one can di- 
rectly predict the outcome of any measurement on the 
system through a simple overlap integral, much the same 
as how we use the Wigner function for predictions. An 
advantage that this method has over other state determi- 
nation techniques is that the amount of state disturbance 
can be minimized. Thus, in principle, we can character- 
ize quantum states in situ, for instance, in the middle of 
quantum computation circuits, or during chemical reac- 
tions, without disturbing the process in which they fea- 
ture. 

We thank Aephraim Steinberg for useful discussions. 
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